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The relativistic approach to electroweak properties of two- 
particle composite systems developed in Ref. [1] is generalized 
here to the case of nonzero spin. In developed technique the 
parametrization of matrix elements of electroweak current op- 
erators in terms of form factors is a realization of the Wigner- 
Eckart theorem on the Poincare group and form factors are 
reduced matrix elements. The p meson charge form factor is 
calculated as an example. 

PACS number(s): 13.40.-f, 11.30.Cp 

A new relativistic approach to electroweak properties 
of composite systems has been proposed in our recent 
paper [1] . The approach is based on the use of the instant 
form (IF) of relativistic Hamiltonian dynamics (RHD). 
The detailed description of RHD can be found in the 
review [2]. Some other references as well as some basic 
equations of RHD approach are given in Ref. [1]. 

Now our aim is to generalize the approach to composite 
systems of two particles of spin 1/2 with nonzero values 
of total angular momentum, total orbital momentum and 
total spin. The main problem is a construction of electro- 
magnetic current operator satisfying standard conditions 
(see, e.g., Refs. [1,3]). 

The basic point of our approach [1] to the construc- 
tion of the electromagnetic current operator is the gen- 
eral method of relativistic invariant parameterization of 
local operator matrix elements proposed as long ago as 
in 1963 by Cheshkov and Shirokov [4]. This canoni- 
cal parametrization of local operators matrix elements 
was generalized to the case of composite systems of free 
particles in Refs. [5,6]. This parametrization is a real- 
ization of the Wigner-Eckart theorem for the Poincare 
group and so it enables one for given matrix element of 
arbitrary tensor dimension to separate the reduced ma- 
trix elements (form factors) that arc invariant under the 
Poincare group. 

Physical approximations that we use in our approach 
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are formulated in terms of reduced matrix elements, for 
example, the well known relativistic impulse approxima- 
tion. In our method this approximation does not violate 
the standard conditions for the current. 

In the present paper we propose a general formalism for 
the operators diagonal in the total angular momentum. 
The details of calculations can be found in [3] . 

Let us consider the operator = j M (0) that describes 
a transition between two states of a composite two- con- 
stituent system. Let us neglect temporarily for simplicity 
the conditions of self-adjointness, conservation law and 
parity conservation. The Wigner-Eckart decomposition 
of the matrix element has the form [4]: 

(Pc mjc \jfi\fc',m' Jc ) = {mj c \D Jc {p c , p' c ) 
x [Ff A; + F^YM) +Fi R» + K K*\ Kc) . (1) 

^ C = E fin(Q 2 )(iPc^{p' c )r . (2) 



n=0 



Here 



£fiv\p P c P, 



A M = (pc - p'Jn = Qf^K'^ = {p c + 
^(p' c );{p c -p' c f = -Q\p\ 



= P'c 2 = 

M c 2 , M c , J c are the mass and spin of the composite parti- 
cle, mj c is spin projection, T p {p' c ) is the spin four-vector 
defined with the use of the Pauli-Lubansky vector [1] , f[ n 
are reduced matrix elements, e^ u \ p is a completely anti- 
symmetric pseudo-tensor in four dimensional space-time 
with e i23 = -L 

In the frame of RHD the form factors of composite 
systems ff n are to be expressed in terms of RHD wave 
functions and constituents form factors. 

In RHD a state of two particle interacting system is 
described by a vector in the direct product of two one- 
particle Hilbert spaces (see, e.g., Ref. [1]). So, the matrix 
element in RHD can be decomposed in the basis [1] 



P, A J, I, s, 



mj , 



(3) 



HereP M = { Pl +p 2 )u,P^ 



s, i/s is the invariant mass of 
the two-particle system, I is the orbital angular momen- 
tum in the center-of-mass frame (cm.), S is the total 
spin in the cm., J is the total angular momentum with 
the projection raj . 
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(Pc mj\j^\p c ', m'j) = £ J ^Jft, d Vs dVs' 

x (p c , m Jc \P, y/s, J, I, S, mj) 
x (P, yfs, J, I, S, mj\j^\P', VP, J', I', S', mj,) 

x(P',^,J',l',S',mj,\p c ',m' Jc ) . (4) 

Here the sum is over variables J ,</' ,1,1' ,S,S' ,mj,mji , and 
(P' ,V~s', J 1 , 1' , S' ,mji\p' ,m' Jc ) is the wave function in 
the sense of IF RHD. 



(P ,\/s,J,l,S,mj\p c ,mj c ] 

N C 8(P - Pc)Sj c jSm Jc mj fl§( k ) 



(5) 



Here k = V /A(s , A/f , M|)/(2 v / s), M x , M 2 are masses 
of constituents, X(a, b, c) = a 2 + b 2 + c 2 — 2(ab+bc + ac). 

The RHD wave function of constituents relative motion 
with fixed total angular momentum is defined as 



(6) 



(7) 



<p{§(k{s)) = \J^(l- V 2 / S 2 )u ls (k)k, 
and is normalized by the condition 

£ f u 2 s (k)k 2 dk = l. 

IS J 

Here 77 = M 2 — Mf , uis(k) is a model wave function. 

The main difficulty arising in this case is the follow- 
ing. In the expression (1) we were dealing with the 
paramctrization of local operator matrix elements in the 
case when the transformations of the state vectors and 
of the operators were defined by one and the same rep- 
resentation of the quantum mechanical Poincare group. 

A different situation takes place in the case of the ma- 
trix element in the r.h.s. of Eq. (4). The operator de- 
scribes the system of two interacting particles and trans- 
forms following the representation with Lorentz boosts 
generators depending on the interaction [1]. The state 
vectors physically describe the system of two free par- 
ticles and present the basis of a representation with 
interaction-independent generators. So, the Wigncr 
Eckart decomposition can not be applied directly to the 
matrix element in the integrand in the r.h.s. of Eq. (4). 
This is caused by the fact that it is impossible to con- 
struct 4-vectors describing the matrix element transfor- 
mation properties under the action of Lorentz boosts 
from the variables entering the state vectors (contrary 
to the case of, e.g., Eq. (1)). In fact, the possibility of 
matrix element representation in the form (1) is based 
on the following fact. Let us act by Lorentz transforma- 
tion on the operator U~ 1 (A)j fJ -U(A) — j^. We obtain 
the following chain of equalities: 



(Pc,m Jc \j tl \p c ',m'j c ) 
= (p c , mjcp-^A^UiAM ', m' Jc ) 
= Y, {rnjoWD^R^-^mjc) 

rhjc ,7ri j 

x (Ap c , m Jc \f\Ap c ', m'j c )(m' Jc \D J *(R A )\ m' Jc ). (8) 

Here D Jc {R\) is rotation matrix realizing the angular 
momentum transformation under the action of Lorentz 
transformations. The equalities (8) show that the trans- 
formation properties of the current as a 4-vector can be 
described using the 4-vectors of the initial and the final 
states. This means that the canonical parameterization 
[4] is the realization of the Wigner-Eckart theorem on 
the Poincare group. 

In the case of the current matrix element in the r.h.s. 
of Eq. (4) the relations (8) are not valid and direct ap- 
plication of the Wigner-Eckart theorem is impossible. 

However, it can be shown that for the matrix element 
in Eq. (4) considered as a generalized function (distribu- 
tion), that is considered as an object having sense only 
under integrals and sums in Eq. (4). So, the equality (8) 
is valid in the weak sense. 

Let us consider the matrix element in question as a 
regular Lorentz covariant generalized function (see, e.g., 
Rcf. [7]). Using Eq. (5) let us rewrite Eq. (4) in the 
following form: 

(Pc,m Jc |j M |p c ',TO' Jc ) 

= E / A^VirfV7^( s )^ c 5 ,( s ') 

l,l',S,S' J 

x(Pc,Vs,J c J,S,m Jc \j^\p c ',Vs',Jc,l',S',m'j c ) . (9) 

Here it is taken into account that the current operator 
is diagonal in total angular momentum of the composite 
system, M = N c N'JN CG N' CG . 

Let us make use of the fact that the set of the states 
(3) is complete: 



dp 

Ncg 



dy/s 



x\P,y/s,J,l,S,mj){P,\/s,J,l,S,mj\ . 



(10) 



Here the sum is over the discrete variables of the basis 
(3). 

Under the integral the matrix element of the trans- 
formed current satisfies the following equalities ((5) and 
(10) are taken into account): 

£ / Nd^d^^( s )ri, c s ,( s ') 
x(p c ,^,J c ,l, S, m Jc \U- 1 {A)j t ,U{A) 
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x\p c ',Vs 7 ,J c ,l',S',m' Jc ) 
= E/ Nd^- S d^rf§{ s ) v j v c s ,{ s ') 

x Y, (mjc\[DHR*)]~ 1 \™Jc) 

Thjc,'rh / j c 

x (Ap c , y/s,J c ,l,S, m Jc \j^\Ap c ', Vs 7 , J c ,l',S', m' Jc ) 

x(m' Jc \D J c(R A )\m[ Jc ) . (11) 

It is easy to see that under the integral the current ma- 
trix element satisfies the equalities analogous to Eq. (8), 
so now it is possible to use the parameterization under 
the integral, that is to use the Wigner-Eckart theorem 
in the weak sense. The r.h.s. of Eq. (9) can be written 
as a functional on the space of test functions of the form 
(see Eq. (6), too)) ^«' ss '( s , s') = u ls (k(s)) Ul < s >(k(s')) , 
and Eq. (9) can be rewritten as a functional in R 2 with 
variables (s,s'): 

{Pc,rn Jc \jn(0)\p c ' ,m' Jc ) 
= E / dKs,s')M^ u ' ss '(s,s') 

l,l',S,S' J 

x(Pc,Vs,JcJ,S,m Jc \j^\p c \Vs',Jc,l',S',rn'j c ) . (12) 

Here the measure is chosen with the account of the rel- 
ativistic density of states, subject to the normalization 
(6), (7): 

dn(s, s') = 16 0(s - (Afi + M 2 f) 6(s' - (Mi + M 2 ) 2 ) 

x y] Vi(l - ry 2 /s 2 )V7(l - V 2 /s< 2 ) dfi{s) dfi{s') . (13) 

Here dfi(s) = (l/4)kdy/s. 

The sums over discrete invariant variables can be trans- 
formed into integrals by introducing the adequate delta- 
functions. The obtained expressions are functional in 
R 6 . 

The functional in the r.h.s. of Eq. (12) defines a 
Lorentz covariant generalized function, generated by the 
current operator matrix element. 

Taking into account Eq. (11) we decompose the matrix 
element in the r.h.s. of Eq. (12) into the set of linearly 
independent scalars entering the r.h.s. of Eq.(l): 

N{Pc Vs, J c , I, S, mjclj^pc ', VP, J c , I', S', m' Jc ) 
= (m Jc \D J °( Pc , p' c )Y(iPc^(Pc)) n 

71 = 

xA l ^ ss '(s,Q 2 ,s')\m' Jc ). (14) 

Here A^^ (s,Q 2 ,s') is a Lorentz covariant generalized 
function. 



Making use of Eq. (14) and comparing the r.h.s. of 
Eq. (1) with Eq. (12) we obtain: 

]T / d^ s , s ')^' ss '( s , s ') 

l,l',S,S' J 

x{m Jc \A l ^ s '{s,Q 2 ,s')\m'j c ) 

= (m Jc \[f^Kl + f^T,(p' c ) 

+fLR l * + fLK l *]\m'jc)- (15) 

All the form factors in the r.h.s. of Eq. (15) are nonzero 
if the generalized function A contains parts that are di- 
agonal (At) and non-diagonal (A2) in m,j c , m' Jc For the 
diagonal part we have from Eq. (15): 

E / dn{s,s')^ U ' SS ' (8, 8') 
l,l'S,S> 

(m Jc \Atf(s,Q 2 ,s')\m Jc ) 
= (m Jc \ [A c n [V] K + / 4 C „[V] K„] \m Jc ) . (16) 

The notation ff n [^] in the r.h.s. emphasizes the fact that 
form factors of composite systems are functional on the 
wave functions of the intrinsic motion and so, on the test 
functions. 

Let the equality (16) be valid for any test function 
ip n ss (s,s'). When the test functions (the intrinsic mo- 
tion wave functions) are changed the vectors in the r.h.s. 
are not changed because according to the essence of the 
parametrization (1) they do not depend on the model for 
the particle intrinsic structure. So, when the test func- 
tions are varied the vector of the r.h.s. of Eq. (16) re- 
mains in the hyperplane defined by the vectors , K^. 

When test functions are varied arbitrarily the vector 
in l.h.s. of Eq. (16) can take, in general, an arbitrary 
direction. So, the requirement of the validity of Eq. (16) 
in the whole space of our test functions is that the l.h.s. 
generalized function have the form: 

Alf'(8,Q 2 ,8') = K' fl Gl ss '(8,Q 2 ,8') 

+K^Gt SS \s,Q 2 ,s') . (17) 

Here G l l n ss (s,Q 2 ,s') , i = 1,4 are Lorentz invariant 
generalized functions. Substituting Eq. (17) in Eq. (16) 
and taking into account Eqs. (6) and (13) we obtain the 
following integral representations: 




xGl ss '( S ,Q 2 ,s>) (18) 

for i = 1, 4. In the case of matrix element in Eq. (15) non- 
diagonal in m Jc , m' Jc we can proceed in an analogous 
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way and obtain an analogous integral representations for 
j-(Q 2 ), * = 2,3. 

So, the matrix element in the r.h.s. of Eq. (12) con- 
sidered as Lorentz covariant generalized function can be 
written as the following decomposition of the type of 
Wigner-Eckart decomposition: 



(p c , y/s, J c , I, S, m Jc \j^\p c ' , Vs 7 , J c , I', S', m' Jc ) 
= ±(m Jc \D J °(p c , p' c ) [Ti K',, + T 2 T»(p' c ) 
-^3 Rfi + ^4 \m' Jc 




0.01 



Q 2 (GeV 2 ) 



FIG. 1. The results of the calculations of the p-meson 
charge form factor with different model wave functions [1,3]. 
The solid line represents the relativistic calculation with the 
wave function of harmonic oscillator, the dashed line - with 
the power-law wave function for n — 3, dash-dot-line - with 
the wave function with linear confinement, dotted line - with 
the power-law wave function for n = 2, dot-dot-dash-line - 
the non-relativistic calculation with the wave function of har- 
monic oscillator. The wave functions parameters are obtained 
from the fitting of p - meson MSR. The sum of quark anoma- 
lous magnetic moments is taken as k u + kj = 0.09 in natural 
units. The quark mass is M — 0.25 GeV. 

In Eqs. (19), and (20) the form factors Gf^ ss ' (s, Q 2 , s') 
contain all the information about the physics of the tran- 
sition described by the operator j^. They are connected 
with the composite particle form factors (1), and (2) 
through Eq. (18). In particular, physical approximations 



are formulated in our approach in terms of form factors 
G^ ss ' (s,Q 2 ,s') (see Ref. [1] for details). The matrix 
element transformation properties are given by the 4- 
vectors in the r.h.s. of Eq. (19). 

It is worth to emphasize that it is necessary to con- 
sider the composite system form factors as the functionals 
generated by the Lorentz invariant generalized functions 

Now let us impose the conditions of self-adjointness, 
conservation law and parity conservation on the matrix 
elements in Eqs. (1), and (19). The r.h.s. of equalities (1) 
and (19) contain the same 4-vectors and the same sets 
of Lorentz scalars (2) and (20), so, to take into account 
the additional conditions it is necessary to redefine these 
4-vectors and functions Gf n ss (s,Q 2 ,s r ). For example, 
the conservation law gives F£ = and J 7 4=0. 

Let us write the parameterization (19), (20) for the 
particular case of composite particle electromagnetic cur- 
rent with quantum numbers J = J' = S = S' = 1, which 
is realized, for example, in the case of deuteron. Separat- 
ing the quadrupole form factor and using Eqs. (19), and 
(20) we obtain the following form: 



(Pc 



, J c , I, S, mj c \jfj,\p c ', V s', J c , I', S', m' Jc ) 



Af 



(mjclD^Pcp'c) 



W i 



L Jc/ 



(21) 



Tx = G%{ s ,Q 2 ,s') + G\\{ s ,Q 2 , s '){[ipc v T v {p'c)\ 2 

T 3 = G%(s,Q 2 ,s / ) . (22) 

We have taken into account that the equation 
G l 2\{s,Q 2 ,s') — is valid in weak sense. The 
parametrization (1), (2) takes the form (21), (22) with 

Gf n (s,Q 2 ,s') —> fin(Q 2 )- It is eas y to see that for 
the redefined form factors the equality (18) remains 
valid. Form factors /f n (Q 2 ) are connected with charge, 
quadrupole and magnetic Sachs form factors: Gc{Q 2 ) — 
koiQ%G Q (Q 2 ) = (2M 2 JQ 2 )K 2 (Q 2 ),Gm(Q 2 ) = 
— M c f 3Q (Q 2 )- The modified impulse approximation 
(MIA) can be formulated in terms of form factors 
G l i q (s, Q 2 , s'). The physical meaning of this approxima- 
tion is considered in detail in Ref. [1,3]. 

The results of calculations for the p - meson charge 
form factor in MIA (I = I' = 0) are represented in Fig.l. 
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problems (see Ref. [8]). 
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